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We study the energy spectra and wavefunctions of graphene rings formed from metallic armchair 
ribbons, near zero energy, to search for properties which may be identified with "effective broken time 
reversal symmetry" (EBTRS). Appropriately chosen corner junctions are shown to impose phase 
shifts in the wavefunctions that at low energies have the same effect as effective flux tubes passing 
near the ribbon surface. Closing the ribbon into a ring captures this flux and yields properties 
that may be understood as signatures of EBTRS. These include a gap in the spectrum around zero 
energy, which can be removed by the application of real magnetic flux through the ring. Spectra 
of five and seven sided rings are also examined, and it is shown these do not have particle-hole 
symmetry, which may also be understood as a consequence of EBTRS, and is connected to the 
curvature induced in the system when such rings are formed. Effects of deviations from the ideal 
geometries on the spectra are also examined. 



PACS numbers: 73.23.-b, 73.23.Ra, 81.05.Uv 



I. INTRODUCTION 

Graphene is one of the most interesting two dimen- 
sional electron systems to become experimentally avail- 
able in the last decade The low energy physics of 
this system is dominated by two Dirac points per spin, 
through which the Fermi energy passes when the system 
is nominally undoped 0. Among the many fascinating 
ideas associated with this system is a concept of "effec- 
tive time- reversal symmetry breaking" Q . There are two 
ways in which this occurs. Most directly, each individ- 
ual Dirac point of graphene does not reflect time-reversal 
symmetry by itself: time reversal in this system involves 
interchanging the two valleys. This can have dramatic 
consequences. For example, transport in graphene sys- 
tems where disorder is ineffective at intervalley scattering 
is expected to lead to weak anii-localization, rather than 
localization @, i, 0, i, 0, fl @ . 

A second related effect occurs when a graphene sur- 
face has curvature. The effective long-wavelength theory 
of the system in such cases include gauge fields which 
can be understood as effective magnetic fields 0, [Io|, UK • 
These fields are directed in opposite directions for each 
valley, so that the system as a whole respects time re- 
versal invariance. If there is no source of short-range 
scattering in the system, such that intervalley scattering 
is negligible, the electrons can behave as if they are mov- 
ing in a non-uniform magnetic field [3[ . Such a viewpoint 
can be used to infer the single particle energy spectrum 
of a buckeyball [10], or in the vicinity of a disclination in 
an otherwise perfect graphene lattice [ill ]. 

Clearly, in the absence of applied magnetic fields 
and/or magnetic impurities, and if interactions may be 
neglected, then the properties of graphene will be time- 
reversal invariant. Nevertheless, one may ask what prop- 



erties can be understood as reflecting effective broken 
time-reversal symmetry in an idealized situation. One 
way to do this is by studying quantum rings. In an ap- 
plied magnetic field, a generic one-dimensional system 
formed into a ring carries a persistent current [l2j . Sys- 
tems with broken time reversal symmetry may also carry 
spontaneous currents in zero magnetic field, as may be 
the case for example in SrRu04 [131 ] . One can thus exam- 
ine the low-energy electronic structure of a graphene ring 
to see what properties it may have in common with sys- 
tems where time reversal symmetry is genuinely broken. 
This is the subject of our study. 

Graphene rings have been studied by several groups. 
Nanoribbons closed into short nanotubes have spectra 
which are sensitive to the precise boundary termination 
of the ring, and may behave in a complicated way as 
the width is varied 0, [l5| ■ Some studies adopt simpli- 
fied boundary conditions allowing for studies using the 
Dirac equation, so that one may construct spectra from 
those of the individual valleys [161 ]. or treat many-body 
effects [17j , Tight-binding studies of flat graphene rings 
involve different edge terminations and corner geometries 
fl6L [HI [T9L [20j . and reveal spectra which are sensitive 
to both, consistent with previous studies of transmission 
through polygons and junctions [2l|, |2^|. Some experi- 
mental studies have recently been reported [23|, [24|, in 
which Aharonov-Bohm oscillations in graphene rings are 
reported. This demonstrates that sufficiently large phase 
coherence lengths may be reached to allow quantum co- 
herence effects to be observed, although current sample 
geometries are not sufficiently controlled to allow direct 
comparison with the theoretical studies of idealized mod- 
els. 

In this work, we will focus on graphene rings con- 
structed from metallic armchair ribbons. The motivation 
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for this choice is that it offers clear contrasts from what 
is expected if the unique properties of its two valleys and 
sublattices are not properly accounted for. Another im- 
portant motivation for this choice is that 60° corner junc- 
tions may be constructed for such ribbons which are per- 
fectly transmitting [2l[ in the lowest transverse subband. 
Under such circumstances one may expect a conventional 
semiconductor to have a state at zero energy, since the 
transverse confinement energy vanishes and backscatter- 
ing at the junctions is absent. For graphene, we shall 
see that while there is no backscattering, there is a phase 
shift upon passing through such junctions which open a 
gap around zero energy even in the absence of disorder. 
This phase shift is highly analogous to one associated 
with the holonomy arising in a "cut and paste" proce- 
dure used to create a disclination; in particular it may 
be modeled by introducing a (singular) gauge potential 
[Tol IfO ] which, in the present case, induces an effective 
flux contained in ring. Thus the repulsion of the spec- 
trum away from zero energy in the absence of any real 
magnetic flux threading the ring may be understood as a 
manifestation of "effective time-reversal symmetry break- 
ing" (ETRSB). 




(a) 




(b) 



FIG. 1: Armchair hexagonal ring, (a) Illustration of a ring 
with inner radius r a , outer radius r&, and flux $ passing 
through the hole, (b) Energy spectrum near E = in units 
of the hopping matrix element t, for E > 0. For this ring 
the spectrum is particle-hole symmetric, and r a = 32.5a, 
rt = 38.5a. Integers label the energy levels. 



Fig. 1(a) illustrates a simple geometry, a hexagonal 



ring constructed of metallic armchair nanoribbons, and 
the associated energy spectrum as a function of mag- 
netic flux through the hole of the ring is presented in 



Fig. 1(b). As we shall see, the spectrum for states closest 
to zero energy may be understood in detail by account- 
ing for the phase jumps at the junctions as well as the 
six-fold symmetry of the system. Levels with different 
rotational quantum numbers cross (without level repul- 
sion) at higher energies in the absence of perturbations 
that allow them to admix. (Similar behavior has been 
observed previously for related ring geometries [lril. [l8|.) 
Studies [16| of zigzag hexagonal rings have noted that 
their energy levels may be interpreted in terms of over- 
laid spectra for each valley, each of which is asymmetric 
with respect to real magnetic flux through the ring, but 
which together yield a spectrum which is symmetric with 
respect to time reversal. Although states of armchair rib- 
bons do not have valley index as a good quantum number, 
we shall see that there is an analogous quantum number 
that can be associated with the states of the lowest sub- 
band, so that in this case too the spectrum may be un- 
derstood as overlaid spectra which individually are not 
time-reversal symmetric. 

Curvature may be introduced by considering rings with 
fewer or greater than six sides. Such systems may be 
related to hexagonal rings by a cut and past procedure 
in which one either cuts out one or more sides of the 
ring, and stitches together the dangling bonds, or cuts 
open the system at a corner and inserts one or more extra 
sides. This is highly analogous to the procedure by which 
one creates a disclination in a perfect graphene sheet [ToL 
[ll[ . Since the procedure results in changing the number 
of corners, the effective flux through the ring changes, 
modifying the spectrum near zero energy in a way which 
may be understood in detail. A prominent feature of the 
spectra of such rings is that they are not particle-hole 
symmetric. In general, one may show that a honeycomb 
lattice in which closed loops always have an even number 
of steps must have a particle-hole symmetric spectrum. 
Interestingly, for five-fold and seven-fold rings, any small 
loop must have an even number of steps, but loops that 
surround the hole of the ring will be odd in number. 
Thus the breaking of particle-hole symmetry is associated 
with a topological property of the ring, and is inextricably 
tied to the curvature induced in forming pentagons and 
heptagons from graphene. Thus the broken particle-hole 
symmetry in the spectrum may be understood as another 
manifestation of EBTRS. 

Clearly, to observe these effects directly requires ex- 
treme control of the lattice, beyond the capabilities of 
current fabrication technologies. We briefly examine the 
effects of disorder and perturbations near the corners and 
edges to see how disturbing these are to the features de- 
scribed above. We find that the most profound changes in 
the spectrum occur when the corners are modified from 
their ideal form: in particular such defects can induce 
zero energy states which are absent in a perfect ring ge- 
ometry. These are localized states which are insensitive 
to the flux through the ring, and so do not induce a 
persistent current. Perturbations near the edges shift 
the spectrum and may contribute nearly localized states, 
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but do not in general bring states to zero energy at zero 
flux, so that the main feature which may be ascribed to 
EBTRS remains. 

This article is organized as follows. In Section [IT] we 
discuss the subband structure of the armchair ribbon, 
and the phase jump associated with creating a 60° corner 
in such a ribbon. Section lTTTl describes our results for ideal 
rings, and we demonstrate that they can be understood 
in great detail from the results of the previous section. 
In Section IIVI we examine the effects of a few deviations 
from the perfect ring geometry, and finally we conclude 
with a summary in Section fVl 

II. METALLIC ARMCHAIR NANORIBBONS 
AND THE TRANSMITTING CORNER 
JUNCTION 




FIG. 2: Geometry for an armchair nanoribbon corner junc- 
tion. Note that atoms from both sublattices appear with 
equal frequency at the edge of each ribbon, implying vanish- 
ing boundary conditions for both. Transverse modes of the 
ribbons may be matched at the junction by equating wave- 
functions on open and closed circles, and on triangles. The 
integers spanning the width of the incoming and outgoing 
ribbons correspond to the index I in Eq. [7] 

I 



where the upper (lower) entry represents the amplitude 
on the A (B) sublattice, W is the ribbon width, K = 
0), K' = (|f , 0), and a is the lattice constant. The 
value of p x must be chosen such that the total amplitude 
at the edges of the ribbons vanishes, so that p x — > p n 
comes in quantized values (25| . For metallic ribbons, the 
lowest subband satisfies p n =o = 0. These wavefunctions 



Figure [2] illustrates two armchair nanoribbons meeting 
at 60° corner junction. A careful comparison of the wave- 
functions and energy spectra of armchair nanoribbons as 
found from numerical tight-binding calculations and from 
analytical solutions of the continuum model based on the 
Dirac equation demonstrates that these may be brought 
into quantitative agreement, provided one adopts van- 
ishing boundary conditions for both the A and B sublat- 
tices at the edges of the ribbon [25|. In the continuum 
this means that the boundary condition admixes states 
from the two valleys around the Dirac points K and K'. 
This is in contrast to the situation for zigzag edges, for 
which the wavefunctions vanish on only one of the two 
sublattices, and there is no valley mixing due to bound- 
ary conditions [25]. Moreover, in the lowest transverse 
subband for a zigzag nanoribbon, the direction of cur- 
rent is associated with a particular valley. This chirality 
underlies proposals to use zigzag graphene nanoribbons 
as elements in "valleytronics" applications [26(. Because 
armchair ribbons admix valleys one may suppose that 
such physics is irrelevant to them. However, armchair 
nanoribbons turn out to support an analogous chirality, 
as we now demonstrate. 



In the continuum limit, the positive energy wavefunc- 
tions of an armchair nanoribbon may be written in the 
form 



(1) 
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have energy s = vp\p\, where vf is the speed of electrons 
near the Dirac points. Negative energy wavefunctions are 
related to those of Eq. [I] by changing 1 — > — 1 in the A 
sublattice amplitudes. 

The metallic case (p n = 0) allows an interesting sim- 
plification of the transverse wavefunction, yielding 
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The subscripts K,K' are introduced in the above ex- 
pression to denote the fact that the spinors provide am- 
plitudes for each of the two valleys. Writing the am- 
plitudes for the wavefunction in vector form, *& Pn ,p y = 
( $ i,K. ]$ s,K* |$ i,K'* ]$ fl,K'.)t i t he p n = wavefunc- 
tion is an eigenstate of the matrix 



T = 



( i 
-i 
i 
V -i 



The matrix T interchanges the amplitudes between 
both sublattice and valley, and multiplies by a phase ±i 
[27I . This particular matrix may be shown to commute 
with the Hamiltonian in the continuum limit, so that 
it represents a symmetry of graphene at low energies. 
Moreover, the boundary condition for armchair edges 
does not violate this symmetry. Thus, any eigenstate of 
a graphene armchair nanoribbon may be expressed as an 
eigenstate of T. In general this requires admixing states 
of different signs of p n . However for the special case of 
p n = 0, one finds 



T^o, Py = sgn{s)sgn{p v )^Q, 1 



Thus, for metallic nanoribbons the eigenvalue of T in 
the lowest subband (p n = 0) is tied to the direction of 
current sgn(p y ), in a way that is highly analogous to the 
connection between current direction and valley index for 
zigzag nanoribbons, with the eigenvalue of T playing the 
role of valley index. 



(3) For our present study of armchair rings, the corner 
junctions as illustrated in Fig. [5]are used. Such junctions 
were shown previously to support perfect transmission 
[2ll ] . The transmission amplitude may be computed by 
matching wavefunctions and currents for the two ribbons 
on the joining surface. Current matching is in general in- 
convenient because it involves products of wavefunctions. 
This requirement may be simplified for the lowest sub- 
band, because the wavefunctions vanish on the open cir- 
cles in Fig. [21 Thus one only need match the currents on 
the bonds connecting the closed circles to the triangles. 
This may be accomplished straightforwardly by match- 
ing the wavefunctions on the triangles as well. For low 
energies, in which we focus on the lowest transverse sub- 
band and take the p y — > limit, referring to the labeling 

(4) in Fig. [2] the resulting matching conditions are [2l| 
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The transmission amplitude between states of the rib- 
bons can be computed in the single mode approximation 
(SMA) 0,H| by computing the overlaps M 0<0 (p y ) of 
the wavefunctions on the joining surface, 



M 0i0 (p y ) = J dx4]l v (x(X),y(X))^ J y (x(X),y(X)), 

(6) 

where ^q 1 ^ ^ s ^ ne wavefunction to the left of the junction, 



(2)*, 



(2) 

■0q is the wavefunction to the right of the junction, 
and A parameterizes the joining surface. Note that in 
the limit p y — > 0, there is no actual y dependence in 

4o( x ->y)- F rom Fig- 121 one may see that the positions 
denoted as x(') = £(') demarcate increments of length 
a/2. The meaning of the formal expression (Eq. 0, 
using Eq. [TJ then takes the form for p y > 0, e > 0, 
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The denominator in Eq. [7] comes from normalizing the 
wavefunctions. Eq. [7] yields Mofi(p y — * 0) = i; more 
generally, one may verify that M 00 (p y — > 0) = z for 
p y e > 0, and M 0j o(p a — > 0) = — i for p y e < 0. The result 
confirms the unit transmission found previously for these 
junctions [21] , and also implies a remarkable property: 
the wavefunctions should jump by a factor of ±i as one 
passes through a junction, with the sign determined by 
the eigenvalue of the T operator defined above. 

This phase shift can also be viewed from the continuum 
perspective. Consider a long metallic armchair nanorib- 
bon segment with its front and back ends identified to 
form a cylinder, effectively a short, fat nanotube. We 
may simulate the phase shift associated with the junction 
shown in Fig. d]by adding a gauge field. Our Hamilto- 
nian has the form (h — 1) 



H( 



VF 
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(8) 

where p x ^ = \d x{y) . Writing * p „, Ph = e lx ^' Pn , we 



can rewrite the Dirac equation for the wavefunctions in 
the form H(p - A)*' = eV' , with A = d y xv 



and x — where the junction between ribbons 

is located at y — 0. In this representation, *f?' pn Py is 
continuous across the junction, and the phase jump is 
fully implemented by the e %x factor. 

The presence of the gauge field can be interpreted as 
being due to a pair of solenoids carrying magnetic flux 
in opposite directions, one above the plane of the ribbon, 
the other below. In this way, one can understand the 
problem of n-sided rings constructed from metallic arm- 
chair ribbons and corner junctions such as shown in Fig. 
[3] as being the same (near zero energy) as the problem of 



® 



FIG. 3: (Color online.) Geometry for flux tubes, arranged to 
produce the effects of the phase factor created by the corner 
junction. Effective flux above and below the plane of the 
ribbon run in opposite directions. 



a ribbon closed into a cylinder (i.e., a short nanotube), 
with In flux tubes, half threaded through in one direction 
and half just outside it in the other direction. In this way 
the system has properties illustrating EBTRS. Just as in 
examples of this effect in graphene [loL [ll[ , real time re- 
versal symmetry is preserved for the system, in this case 
because the effective flux runs in opposite directions for 
different eigenvalues of the operator T. As we shall see 
in detail below, the low energy spectra of such rings as 
found from solutions of the tight-binding model behave 
precisely as if these phase jumps are present. The result- 
ing spectra present properties which may be understood 
as signaling the EBTRS in graphene. 



III. NUMERICAL RESULTS 

The simplest ring system one can study using metallic 
armchair ribbons and the 60° corner junctions discussed 



above is the hexagonal ring, as illustrated in Fig. 1(a) 
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The fact that the junctions are perfectly transmitting in 
the lowest subband might lead one to think that the low 
energy spectrum is the same as that of a metallic arm- 
chair ribbon closed into an cylinder (i.e., a very short 
carbon nanotube.) If this were the case, one would ex- 
pect states at zero energy when no external magnetic 
flux threads the ring. Our discussion above however in- 
dicates that one needs to include the effective flux passing 
through the ring to understand the spectrum. 

Fig. 1(b) is the spectrum obtained from computing 



the eigenvalues of the tight-binding model near zero en- 
ergy, as a function of flux $ through the ring. Note in 
these calculations we include only the phase factors in the 
hopping matrix elements due a solenoid passing through 
the hole of ring; magnetic flux through the individual 
plaquettes of the honeycomb lattice is not included. In 
general the spectrum in this case is perfectly particle-hole 
symmetric, so only positive energy states are displayed. 

The form of the spectrum is highly reminiscent of what 
has been seen previously in hexagonal rings with zigzag 
edges 16]. In particular it takes the form of two sets of 
spectra, each with broken time reversal symmetry such 
that the energies have a particular sign of slope near 
$ = 0. The spectra are effectively time-reverses of one 
another so that the spectrum as a whole has time-reversal 
symmetry; in particular the spectrum evolves in the same 
way whether positive or negative flux is threaded through 
the hole. In the case of zigzag ribbons the two sets 
of spectra are associated with the two valleys. In the 
present case, they are associated with the two eigenval- 
ues of the matrix T. Note that the crossing of the energy 
states through zero at $ = ±<&o/2 may be understood as 
resulting from the sum of the effective fluxes due to the 
corner junctions and that of the real field summing to an 
integral number of flux quanta. 

The phase jumps associated with the corners may be 
demonstrated explicitly by a careful examination of wave- 
functions. Fig. 2] labels a set of sites around the ring, 
and Fig. O illustrates the wavefunctions on one of the 
sublattices for the four lowest positive energy levels at 
$ = 0. The jumps in amplitude associated with pass- 
ing through the corners of the junctions are quite appar- 
ent. Moreover, the amplitudes in the sides of the ring 
reflect the e %VvV phase factor in Eq. Q] very accurately, 
provided one recognizes that diagonalization of the tight- 
binding Hamiltonian results in real eigenvectors, so that 
our wavefunctions are actually linear combinations of 
states with positive and negative values of p y . To see 
this, we fit the wavefunctions to functions of the form 
ip(y) = Asm[p y y + <fio], where A and 4>o are constants, 
but inserting phase jumps of ±f at the discontinuities in 
the wavefunctions. Typical results are illustrated in Fig. 
[5J One may see the excellent fit to wavefunctions with 
well-defined p y . We have confirmed that similar phase 
jumps and fitted wavefunctions are produced at least to 
the ninth energy level. 

The results also directly support the expected sign of 
the phase jump, associated with the eigenvalue of T. To 




FIG. 4: Hexagonal ring illustrating site labels for examina- 
tion of wavefunction. Black dots indicate sites on A sublat- 
tice, gray dots are on B sublattice. The actual ring used in 
calculations has r a = 32.5a and = 38.5a, with a the lattice 
constant of the underlying triangular lattice. 
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FIG. 5: (Color online.) Wavefunctions at zero flux $ for the 
four lowest eigenstates n = 1, 2, 3, 4 illustrated in Fig. 1(b) 



see this, we recognize that the wavefunctions are linear 
combinations of plane waves in y of equal and opposite 
p y , such that one obtains real eigenvectors. The val- 
ues of p y are determined by requiring for the hexagon 
that e (w y L+iTr/2)x6 _ ^ wncrc L is the length of one of 
the hexagon sides, and p y > 0. This requirement very 
accurately produces the values of p y used in fitting to 
sine functions, such as in Fig. El as well as the energies 
e = vp\p y \ (see Table I.) The actual wavefunctions are 
linear combinations of e ± (' lPyy+l7TS / 2 \ with s an integer 
s = 0, .., 5 labeling which leg of the hexagon y is on; thus 
we see the phase jump must have different signs for the 
two signs of p y . Finally, we have examined the states 
with e < 0, and indeed find a change in sign for the 7r/2 
phase jumps relative to wavefunctions with the same p y 
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Energy level 


PyL 


0o 


m 


1 


vr/6 


tt/2 
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-7r/6 


-tt/2 


-2 


3 


tt/2 


tt/2 
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-tt/2 


-tt/2 


-3 
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5tt/6 


tt/2 


4 


6 


-571-/6 


-tt/2 


-4 



TABLE I: Table showing values of p y , 9o, and m in formula of 
form 6(p y L + 8q) — 2nm used in matching numerically gen- 
erated wavefunctions with forms expected from wavefunction 
continuity around a ring (see text.) 



and £ > 0. 




(a)n=l; A sublattice (b)n=3; A sublattice 



FIG. 6: (Color online.) Wavefunctions at zero flux $ for the 
n — 1,3 eigenstates with phase jumps added (see text) and 
fitted to sine functions. 

The presence of six tt/2 phase jumps (due to the six 
corners of the hexagon) implies that p y — is not an 
allowed momentum for the electron wavefunction in this 
type of ring. Thus there is no allowed zero energy state, 
as would be expected for a metallic ribbon closed into an 
annulus, and its absence is imposed by the presence of 
phase factors that are very suggestive of effective mag- 
netic flux penetrating the ring, as discussed above. This 
the gap in the spectrum around e = may be interpreted 
as one signature of EBTRS. 

Rings with numbers of sides differing from six may 
also be considered, such as those illustrated in Fig. [7] 
Because the number of corners in these structures is dif- 
ferent than that of hexagons, the effective flux through 
such rings will be different, leading to extra factors of 
±i which must be included in determining the wavefunc- 
tions. These may be examined by solving the correspond- 
ing tight- binding models, and one finds that the spectra 
and wavefunctions may be fully understood in the same 
way as described above for hexagonal rings. Examples of 
spectra are illustrated in Fig. [51 At $ = 0, the energy 
levels precisely obey e = ±vp\p y \, with p y chosen to be 
consistent with the total phase jump associated with the 
number of corners. Fits analogous to those shown in Fig. 
[5] may be made with equal success. 

The most prominent difference between the spectra 
displayed in Fig. [SJand those of hexagonal rings is the ob- 
vious lack of particle-hole symmetry. In the continuum 
limit, particle-hole conjugates are related by changing 
the relative sign of the amplitudes on different sublat- 
tices (equivalent to multiplying the wavefunction by the 
Pauli matrix cr z for each valley), and the resulting wave- 




(a)Pentagon 




(b)Heptagon 



FIG. 7: Examples of (a) five-sided and (b) seven-sided 
graphene rings with 60° corners, along with points used for 
analyzing wavefunctions. Note such rings necessarily have 
curvature if their bond lengths are undistorted. 



function has the same value of p y , but an energy of op- 
posite sign relative to that of the original wavefunction. 
For geometries containing the junctions studied in this 
work, however, there is a further effect: the phase jump 
for a wavefunction with a given eigenvalue of T changes 
sign (Eq. U) For even numbers of corners this has no 
net effect (e.g., (±i) 6 = —1), but for odd numbers this 
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FIG. 8: Energy spectra as a function of flux for (a) pentagonal 
ring, and (b) heptagonal ring. Ribbons used in constructing 
the polygon edges are identical to those used for hexagonal 
rings with r a = 32.5a and rj, = 38.5a. 



sheets with simple nearest neighbor hopping of uni- 
form magnitude will support particle-hole symmetry if 
all paths beginning and ending at the same point have 
an even number of steps [2!|. Locally this is the case 
for all the structures considered in this paper. However, 
for rings with an odd number of sides, paths that enclose 
the hole of the ring will have an odd number of steps, 
leading to the breaking of particle-hole symmetry. The 
fact that one must examine paths that enclose the hole 
to uncover the symmetry breaking is consistent with the 
idea that the effect is related to flux through the hole in 
the gauge-field description. 

IV. DISTORTED GEOMETRIES 

In this work we have focused on ideal graphene rings 
in order to focus on signatures of EBTRS in their spec- 
tra. However, currently available fabrication techniques 
do not allow the creation of such perfect rings, and an ob- 
vious question is how robust the features discussed above 
might be with respect to deviations from the ideal struc- 
ture. In this section we examine this question in a few 
simple cases. 



affects the wavefunction continuity condition such that 
p y in general cannot take the same values for positive 
and negative energies. In terms of effective flux through 
the ring, this is a reflection of the fact that the associated 
vector potential is a matrix rather than a scalar quantity, 
allowing the system to have different values for the effec- 
tive flux for positive and negative energies. Moreover, 
for a given eigenvalue of T, the effective flux through 
a six-sided ring is effectively ±$o/2, while for five- or 
seven-sided rings it is ±0o/4. This effective flux can be 
canceled by adding real magnetic flux through the ring, 
bringing some states to zero energy, but the amount of 
flux required depends on the number of sides [cf. Figs. 
T^b)|and E] 

Thus, the absence of particle- hole symmetry in the 
spectrum may be understood as a unique effect of the 
types of gauge fields which can occur in graphene. As 
in the case of disclination defects [HI EH , this effect is 
intimately connected to the presence of curvature in the 
system: one may create pentagons and heptagons from 
hexagonal rings with cut-and-paste procedures that leave 
the local 3- fold coordination of the bonding intact, but in 
doing so curvature is necessarily introduced, and an inte- 
gral number of corners must be added or subtracted from 
the system. When this number is odd, particle-hole sym- 
metry breaking will necessarily be present in the spec- 
trum. Thus, this symmetry breaking is intimately con- 
nected to the introduction of curvature associated with 
this type of topological defect. 

Finally, it is interesting to consider how particle-hole 
symmetry may be broken from the perspective of the 
tight-binding model. It is easy to show that graphene 




(a) 




(b) 



FIG. 9: Model of a hexagonal ring with passivated edges and 
associated energy spectrum, (a) Bold bonds have hopping 
matrix elements given by 0.9t. (b) Energy spectrum as func- 
tion of $. r a = 32.5a, r;, = 38.5a. 

The graphene ribbon edges considered in idealized 
models such as considered here have dangling bonds since 
atoms at the very edge only have two nearest neighbors. 
In practice these bonds are often passivated with hydro- 
gen or oxygen, and moreover relax to slightly different 
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bond lengths than in the bulk, so that the gapless spec- 
trum of ideal metallic ribbons develops a small gap [301 ] . 
This can be modeled [3(| by changing the hopping matrix 
element along the edge in the very last row, as illustrated 
in Fig. [9] One may see that the major effect of this is 
to open a large gap around energy E — 0.04t, where for 
the ideal ribbon this is only a small anticrossing. Inter- 
estingly, the gap to the first state at $ = is only very 
slightly suppressed, so that this basic effect is unchanged 
by the edge modification. Moreover, at $ = ^o/^, the 
levels still cross zero energy, behaving exactly as if the 
real flux cancels the effective flux associated with the cor- 
ners. Thus we find that the qualitative effects of EBTRS 
are robust against this edge modification. 
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(b) 

FIG. 10: (Color online.) Hexagonal ring with a single bond 
removed at the edge, r a = 5.5a, rj, = 11.5a. (a) Illustration 
of a defective bonding structure for small size system. Bold 
bond and associated sites on top edge have been removed, (b) 
Energy spectrum as function of <&. 

More disruptive effects occur when the edges or cor- 
ners become disordered. Figure [TU] illustrates the effect 
of removing a single pair of sites on one edge of the sys- 
tem. One may see in the spectrum that all the crossings 
present for the perfect system now become anticrossings. 
This is easily understood: in the case of the ideal system 
there is perfect six-fold symmetry, and each eigenstate 
carries a rotational quantum number m. The perfect 
symmetry implies that when states with different values 
of m approach one another as a function of they can- 
not admix, and the states will cross in energy. Once the 
symmetry has been broken, all states can are generically 
admixed, and the gap openings seen in Fig. [TUJ follow 
naturally. In spite of these gaps, we see the basic struc- 
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FIG. 11: (Color online.) Hexagonal ring with a row removed 
at the edge, r a = 5.5a, — 11.5a. (a) Illustration of a defec- 
tive bonding structure for small size system. Bold bonds and 
associated sites of top edge have been removed, (b) Energy 
spectrum as function of 

ture of the large gap around zero energy which closes as 
$ approaches $o/2 is maintained. 

If a whole row of sites is removed, as in Fig. [TT] the 
spectrum loses its essential features, and a series of states 
near zero energy that vary little with $ are apparent. 
This indicates that the states have been localized; since 
they do not surround the hole of the ring they are insensi- 
tive to the flux through it. This effect may be understood 
if one notes that the upper armchair ribbon is one row 
narrower than the other ribbons, and hence is not con- 
ducting: such a ribbon in isolation has a substantial gap 
near zero energy (2o| . Thus we do not expect the lowest 
energy states to penetrate into the upper arm of the ring, 
and the various effects due to any effective flux quanta 
associated with the corners will be suppressed. 

Perhaps the most surprising effect occurs when a single 
site is removed from a corner, as illustrated in Fig. [12] 
In this case one finds a non-dispersive state precisely at 
zero energy. This localized state appears to be associated 
with the three point zigzag edge at the corner that results 
from removing this site; analogous localized states ap- 
pear in other defect structures as well when they include 
such three point geometries [3l[ . In this case the low en- 
ergy structure associated with the other corner junctions 
is nearly absent, suggesting the defective corner junction 
has become largely reflective. This demonstrates that the 
geometry of the corner junctions plays a crucial role in 
the results associated with EBTRS in graphene rings: a 
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FIG. 12: (Color online.) Hexagonal ring with a site removed 
at one corner, r a — 5.5a, r% = 11.5a. (a) Illustration of a 
defective bonding structure for small size system. Bold bonds 
and associated sites have been removed, (b) Energy spectrum 
as function of $. 



defective junction will generically admix states with dif- 
ferent T eigenvalues, and moreover can effectively block 
transmission at low energies so that effects associated 
with the non-compactness of the geometry are difficult 
to see or absent at low energies. The sensitivity of trans- 
mission to the precise geometry at a corner junction can 
also be seen in the conductance through such junctions 
(2lj . Similar sensitivity to precise geometries has been 
observed in conductance through graphene polygons and 
quantum dots [2l|, [22j • 

We conclude this section with two remarks. One way in 
which the energy spectrum of a quantum ring is probed 
is through the persistent currents it supports in a real 
magnetic field, J = — c4^. Such currents create magne- 
tization which may be detected either for a large enough 
ring or for a collection of identical rings. Since localized 
states have zero or very small |^|, their presence will 
not be apparent in experiments that probe the system 
in this way. However they should be visible in tunnel- 
ing experiments. Secondly, in experiments which probe 
the spectrum through the magnetization of a collection 
of rings, inevitably there will be variations in the precise 
structure of each ring, so that few if any will have the 
ideal spectra of the last section. However one may spec- 
ulate that some of the features may be reflected in the 
total magnetization as a function of Fermi energy, since 
this involves an average over many different disorder re- 
alizations around the ideal ring structure. We leave this 



possibility for future research. 

V. SUMMARY 

In this paper we have studied the states of graphene 
rings formed from metallic armchair ribbons, with 60° 
corner junctions which are perfectly transmitting at low 
energies. We showed that states of the ribbons have a 
quantum number that is the eigenvalue of a symmetry 
operator T, and states with the same value of this near 
zero energy are connected with magnitude one through 
these particular corner junctions. However the ampli- 
tudes support a phase jump as one passes through the 
junctions, and these phase jumps can be described in 
terms of a gauge field due to effective fluxes passing just 
above and below a ribbon. When such a ribbon is closed 
into a ring, there is an effective flux contained such that 
p y , the momentum along the ribbon direction, cannot 
vanish if the wavefunction is to meet appropriate bound- 
ary conditions. As a result there are no states of a six- 
sided ring at zero energy in the absence of an extermal 
magnetic field. States can be driven to zero energy by 
application of such a field, which effectively cancels the 
flux associated with the junctions. This behavior may be 
understood as a manifestation of effective broken time 
reversal symmetry (EBTRS). 

Five-fold and seven-fold rings enclose different total 
effective fluxes, shifting the spectra such that particle- 
hole symmetry is broken. This behavior is another signal 
of EBTRS, and moreover is intimately connected with 
the curvature induced in the rings when these systems 
are formed. Analogous descriptions have been applied 
to graphene sheets with disclinations [13, [HI > although 
the breaking of particle-hole symmetry in the spectrum 
is not so apparent in those cases as in the present one. 

Finally, we examined the effects of several types of vari- 
ations from the ideal ring structure on the spectra. A 
change near the edges to model relaxation and passiva- 
tion was found to only have a small quantitative effect at 
the lowest energies. Removal of a small number of atoms 
from one edge opened small gaps at points where anti- 
crossing occurs in the ideal case, but removing a whole 
row rendered one side of the ring non-conducting, effec- 
tively removing the signature of EBTRS. Removing a 
single atom from a corner introduced a zero energy state 
and greatly reduced the sensitivity of the system near 
zero energy to the flux through the ring. We conclude 
from this that well-formed corner junctions of the type 
described in this work are a necessary ingredient to see 
effects attributable to EBTRS. 
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